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We use the density matrix formalism to analyze the interaction of interferometer-type supercon- 
ducting qubits with a high quality tank circuit, which frequency is well below the gap frequency of 
a qubit. We start with the ground state characterization of the superconducting flux and charge 
qubits. Then, by making use of a dressed state approach we describe the qubits' spectroscopy when 
the qubit is irradiated by a microwave field which is tuned to the gap frequency. The last section of 
the paper is devoted to continuous monitoring of qubit states by using a DC SQUID in the inductive 
mode. 
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I. INTRODUCTION 

The existence of a superposition between macroscopi- 
cally distinct quantum states has important implications, 
in particular, for the quantum measurement theory, since 
a single macroscopic quantum system can be easily ac- 
cessed by the macroscopic measuring device. One of the 
most advanced solid state quantum systems is a super- 
conducting qubit which is based on either the charge 
or phase degrees of frccdomii^. Several ways of reading 
out the qubit properties have been proposed and imple- 
mented. In an initial research stage switching current 
measurements combined with the excitation of a qubit 
by microwave pulses were used^i^. Recently it turns out 
that in some cases inductive (dispersive) measurements 
can be more effective. In the frame of this approach 
the qubit is coupled with linear— '^'^'^ i^^'^^i^^'^'^i^"^ or non- 
linear oscillators^'^ ■■^'^ijJ. 

In this paper we study two kinds of interferometer type 
superconducting qubits (flux and charge qubits) which 
are realized in practice as Josephson junctions embedded 
in a superconducting (interferometric) loop^ i^°'^^ . 

A distinct feature of our approach is that a qubit is 
inductively coupled with a high quality tank circuit (a 
linear LtCt oscillator, where Lt is the inductance and 
Ct is the capacitance of the tank) which resonance fre- 
quency is well below the gap frequency of the qubit. The 
essential information about the qubit properties can be 
extracted from the voltage across the tank. The method 
is found to be a reliable technique for the investigation of 
the properties of the flux qubit systemsiSiiSiiSSiHi, (prior 
to 2004 see the review paper^- and references therein). 

Our approach here is based on the rigorous quantum 
mechanical formulation of the problem, which allows us 
to apply it to the ground state characterization not only 
of the flux qubifeS, but of the charge qubit as well. In 
addition, we show in Section 3 that our method is ex- 
tremely useful for the investigation of the spectroscopic 
properties of superconducting qubits. 

Since the qubit characteristic frequencies are well 



above the tank resonance we may consider the tank cir- 
cuit as a classical system while the qubit is treated quan- 
tum mechanically. This description has been proposed in 
Refs. 0, and the present work is its natural generaliza- 
tion and continuation^^. 

Therefore, the voltage, V , across the tank obeys the 
equation: 



V + 7tV + uj\V 



dt 



u'^LTh, (1) 



where 77^ — ujt/Qt, Qt is the quality factor of the un- 
loaded tank. A/ is the mutual inductance between the 
tank and the qubit loop, ujt = {LtCt)~^^^ is the tank 
resonance frequency, and Ii, (t) = Iq cos ujt is the tank bias 
current with bias frequency to. Here (Iq) is the quantum 
mechanical average of the quantum operator of the cur- 
rent Iq in the interferometer loop. 
Eq. Ill) in Fourier components reads: 



V{lu) {u 



+ iui^r) — —Mijj> 




iuJLOrpLTlo- 

(2) 



In what follows we assume that the interaction between 
the qubit and the tank is weak. In this case it would be 
reasonable to assume (below we prove this assumption 
by detailed calculations) that 




MIt{lo)Z{uj), 



(3) 



where Z{lu) is a complex function {Z{uj) ~ Zi{uj) + 
iZiiio)) which depends on the qubit properties only. As 
is seen from ([3]) the quantity d{Iq)/dt is different from 
zero only due to the interaction of the qubit with the 
tank. 

By taking into account that It{^) = —iV{uj)/u!LT, we 
obtain for the tank detuning ^, and friction r^: 

^ = ^2 _^2_^fc!Viz2(,^), (4) 



2 



k LqLO, 



(5) 



where Lq is the inductance of the qubit loop, fc^ — 
M'^ /LxLq is a couphng constant, which characterizes the 
inductive interaction between the qubit and the tank. 

From ([4]) and ([5]) we obtain the voltage amplitude 
Vt, and the phase a: Vt — ujuj^LtIq/ ^/^'^ + ^^^t^ 
tan a — ^/ujTt- 

Below we calculate in a rigorous quantum mechanical 
way the low frequency linear response of a qubit to a low 
frequency probe signal. It is necessary to say that our 
approach is valid if the interaction between the qubit and 
the tank is weak which allows us to neglect all nonlinear 
terms caused by the finite amplitude of the low- frequency 
signal. Therefore, the main aim here is to calculate the 
quantity Z{uj) which transfers the qubit properties to 
the low frequency characteristics of the tank circuit: its 
detuning (g]) and friction 

The paper is organized as follows. In Section 2 the 
general approach to the ground state characterization of 
qubits is considered. At the end of the section we ob- 
tain the expressions for detuning and friction both for 
the flux and charge qubits. In Section 3 we consider an 
interferometer type qubit which is inductively coupled to 
a low frequency resonant circuit and subjected simulta- 
neously to microwave radiation. By making use of the 
dressed state approach^! we calculate the low frequency 
susceptibilities for charge and flux qubits. The results ob- 
tained in this section are applied to explain the recently 
found phenomenon of the damping and amplification of 
the output signal by a superconducting flux qubitSS. The 
method developed in this paper is applied in Section 4 
to a continuous readout of a superconducting qubit by a 
DC SQUID in the inductive mode. 



II. GROUND STATE CHARACTERIZATION OF 
QUBITS 

In this section we consider a qubit which is isolated 
from microwave fields and interacts inductively with a 
high quality tank circuit. In what follows we assume the 
tank frequency is well below the qubit gap frequency AE 
and, in addition, the temperature is sufficiently low, so 
that ksT << AE. Under these conditions the signal de- 
tected at the output of the tank circuit is mainly defined 
by the properties of the qubit ground state. 



A. Rate equation for two level system 

In the basis of eigenstates |1), |2) the Hamiltonian of 
a qubit reads: 



Ho = —'^z, 



(6) 



where AE is the gap between two energy states. The 
eigenstates are denoted in the following as |1) (ground 
state) and |2) (excited state) with the properties: 
az\l} - az\2) - |2), ax\l) = |2), ax\2) = 

aY\l)^i\2), f7y|2) = 

Next we introduce the density matrix: 



da ^,rT 1 



(7) 



and define its matrix elements as pn = (l|fT|l), P22 = 
(2|ct|2), pi2 = (l|fT|2), P21 = (2|a|l) = p+ . From we 
find the equations for the elements of the density matrix: 



dpi 



dt 

dp22 

dt 



0, 



dpi 



dt 



'iflpi 



(8) 
(9) 
(10) 



where il = AE/h. In the case of damping the above 
equations should be changed to: 



dpi 



dt 

dp22 

dt 

dpi2 
dt 



= -Fipii -f r|P22, 

= Tipii - T1P22, 
= -iQpi2 - T^pi2, 



(11) 
(12) 

(13) 



where Fj is the transition rate from state |2) to state 
|1) (relaxation rate), Fj is the transition rate from state 
|1) to state 1 2) (excitation rate), and the quantity F;^ is 
the rate of decoherence. For equilibrium conditions the 
relaxation and excitation rates are related by the detailed 
balance law: 



Ft ^ F, exp 



From we obtain 



F_ri 



AE 



tanh 



AE 
2kBT 



(14) 



(15) 



where we define F_ = F| — F| and the longitudinal re- 
laxation time T^^ = F| + F|. 

We rewrite Eqs. (fTTI) . IT^ . (|13p in operator form 



(16) 



where 



L=-FT|l>(l|a|l)(l|+Fjl)(2H2)(l| 

+ FT|2)(l|a|l)(2|-Fj2)(2|a|2)(2| 

-F^|l)(l|a|2)(2|-F^|2)(2H1)(1|. (17) 
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As it is seen from pT|) and (|12p . the total population is 
constant: ^(pii + P22) = 0. We take the normalization 
condition as pn + P22 = 1- The rate equations can be 
further simplified by introducing new variables: p — pw — 
P22^ which is the difference in populations between the 
lower and the higher levels, and /9+ = p\2 + /O21, P- = 
P12 - P21- 



where H is given in ((2T|) . 

From ([25|) we get the following equations for the ele- 
ments of the density matrix: 



dp 2iXIx *o 2A/y $0 



dt 



27r 



P- 



h 2n' 



1 1 

Ti Ti d(j)x 

(26) 



dp 1 ^ 



(18) _^,f7p+-r,p_ + ^-p^ 



2iA/x $0 
/i 27r 



P, (27) 



(ip_ 



(19) 



ofp+ 



iiip^ - r^pp^ 



-P- + -^TtP- (28) 



h 27r 



27r 



dp+ 
dt 



-iilp^ — r^p^ 



(20) 



The quantity p*^'^''^ in (fT5)l is just a steady state solution 
of the equation psp . and is the difference in equilibrium 
populations between the lower and the higher levels. 



B. Interaction between qubit and tank circuit 

The Hamiltonian of a qubit which is coupled to the 
tank in the eigenbasis reads: 



(21) 



where It is the current in the tank inductance, $x is a 
DC bias flux through a qubit loop. 

The current operator of the qubit, Iq, can generally be 
written as: 



Iq — Ixf^x + lyfyy + Iz<yz, 



(22) 



C. Linear susceptibilities for the qubit 

We find the solution to Eqs. (pS)) - ((77)) . by assum- 
ing the coupling A is small. Therefore, the time depen- 
dent solution for these equations can be obtained by the 
perturbation method as the small time dependent cor- 
rections to the steady state values: p{t) = p(°) -I- p(^^(t), 

p+{t) = pf +pi')(t) , p-(t) = pL°^ +pW(0, where pW, 

p^\ and p*^-* are steady state solutions for Eqs. ((26)) - 
(l27l) with A 



0. For these steady state values we readily 

obtain: p(o) = p^'"}\ p^°^ = 0, pL°^ = 0. Therefore, the 

fist order corrections p^^^{t), p'^^\t), p^l\t) satisfy the 
following equations (below we omit superscript 1): 



dp _ 1 ^ ^ 1 dp^'^?) 
dt Ti ^ Ti d(j)x 



dp^ 

—j^ = tilp+ - r^p- 



2iXIx $1 
h 2^ 



Jeq) 
H 1 



(29) 



(30) 



where the quantities Ix, Iy, Iz, which will be specified 
below, depend on the nature of the qubit. 

The interaction with the tank also influences the F 
rates in Eqs. (fTT|) . Therefore, we may write in 

linear approximation 



X 



Ti +A- 



dF 



(23) 



(24) 



where A = 2t: M It /^o^ (px = 27r$x/*o, $0 = h/2e is 
the flux quantum. 

As for Fi^ , we assume this rate is sufficiently high which 
allows us to neglect its modification by slow time depen- 
dent external flux MIt- 

The equation for the density matrix is similar to (|16p 



d<7 



(25) 



dp+ 
dt 



= ir2p_ — T^pp+ 



h 2t:' ■ 



(31) 



From these equations it is not difficult to find the linear 
susceptibilities of the system (Xp(w) = p(a;)/A(w), etc., 
where \{lo) — 2ttMIt{i^) /^o)- Therefore, we get 



Xp(^) 



1 



(iwTi + 1) 



2'Khd{uo) 



Xp+ (t^) = 



2$op' 



2Tihd{ui] 



[{iuj + T^)Iy - nix], 



where 



d(w) = {iu + T^Y + 



(32) 



27<I>nfl('^'^) 

Xp-{^) = ^r^7r^[{^^ + r^)Ix+niy], (33) 



(34) 



(35) 
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D. Calculation of Zi{lj) and Z2{uj) 

First we calculate the average current (Iq). By using 
22ll we obtain: 



(/,) = IxTr{aax) + lyTriaaY) + IzTr{aaz), (36) 

where Tr{aaz) ~ — p, Tr{aax) ~ p+, Tr{aaY) = ip-- 
Therefore, 

{Iq) = -Izp + Ixp+ + Hyp-- (37) 
With the aid of equations and ^ we obtain: 

d{Iq) IZ f ^ .V-^')' 

= — [ P + r_ri - A 



dt Ti \ acpx 

+ I (Ixn - IyT^) P- - [IxT^ + lyn) p+. (38) 

If the interaction between the qubit and the tank is ab- 
sent (A = 0), then, p(o) = p^"^ = 0, pL"^ = 0, and, 
as is seen from (I57|) and ([55)1 , the average current is pro- 
portional to the difference of equilibrium populations be- 
tween the two states: (iq) = — /^p''^'^^ but d{Iq)/dt — 0. 
However, if the interaction is sufficiently weak, such that 
A ^ 1, the first approximation for the latter quantity 
gives: 

dt Ti \ d(j)x 

+ i [Ixn - IyT^) pL'^ - [IxT^ + Iy^) p'^l\ (39) 

and for its Fourier component: 



d{Iq) \ h , , . 



dt 



Ti 



dp'-'"''>\ 27rM 

$0 



'>x 



Jt(w) 



+ ^^t(^) (Ixn - IyV^) Xp_ (w) 

- (/xr^+/yf^)Xp+(^)] ■ (40) 
Therefore, we obtain for Z\{ijj) and Zi{ix}) defined in ([3]): 



$0^1 



(7x0 - IyV^) x;'_ {oj)-— {IxT^ -f Iy^) x'p^ {u 



(41) 



^2(^) = —T^rM 



-— (/^f] ~ IyT^) x;_ (^)-— (/xr^ + Iy^) (^), 



(42) 



where x'p{^), x'p_ (^), Xp+ (w) and Xp(w), x'p_ (t^), Xp+ ('^) 
are, respectively, real and imaginary parts of the corre- 
sponding susceptibilities (|32|l - p4|) . 

If we account for the fact that the gap frequency is 
large, SI ^ ^jTi^, we can simplify the calculations of 
Zi{uj) and (^) to obtain the tank detuning ^, and fric- 
tion, Ft in the following form: 



2 2 T 

UJj, — u — Lz 



(1 



d<^x 

,2 



2k LqUj^ (g^-, 
AE ^ 



(43) 



k^UJ^LqTl dp^'"i^ 

(l+a;2T2) 

+ + /!-)• (44) 

The expressions (j43p . (j44p are applicable for any kind of 
interferometer type superconducting qubit once the com- 
ponents Jx, Iy , Iz of the current operator are known. 

The terms in right hand sides of (|43p and (|44p which 
are proportional to dp^^^^^ /d^x reflect the effect of thcr- 
malization which value depends on the interplay between 
the relaxation rate T^^ and the tank frequency uj. 

The quantity dp('=«) can be expressed as: 



dp^'^i') 
d<^x 



cosh 



\2kBT ) 



where we used the fact that Iz can be written as the 
derivative of the ground state energy Eq over the mag- 
netic flux: Iz = dEo/d^x, where Eq = -AE/2. 
Since for a proper qubit operation the condition AE >> 
ksT is necessary, the quantity dp^'^'^^/d^x scales as 
exp ^— -^^y Therefore, no matter what the value luTi 
is, the contribution to the tank response of the terms in 
(|43)) and ((44l) which are proportional to dp^'^'^^ /d^x can 
be neglected as compared with the contribution of the 
last terms in these expressions. 



E. Low frequency response of the flux qubit 



For the flux qubit the energy gap is 



AE = A, = VA2 



(45) 



where A/2 is the tunnelling amplitude between degener- 
ate flux states. The bias e is controlled by an external dc 
flux $x- £ = 2^oIqfx, where Iq is the critical current of 
the flux qubit, fx - ($x/$o - 1/2). 

The current operator of the flux qubit in the eigenstate 
basis is: 



II 
A, 



{eaz - Aax) ■ 



(46) 
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Therefore, the components of the current operator for 
a flux qubit are as follows: 

Ix = -/?A/Ae, Iz = Iqe/A„lY = 0. (47) 

From (|43|) and (|44p we obtain the following expressions 
for the frequency detuning and the friction of the flux 
qubit 

^^cu^^cj' y ' , (48) 



7t + 



I a7 



A3 



r o^"') 



(49) 



where we neglected the terms which are proportional to 

The last term in right hand side of ((48l) , which gives the 
main contribution of the flux qubit to the tank detuning, 
coincides with the result obtained earlier theoreticallj*^ 
and confirmed by experiment^. 

The last term in can be expressed in terms of the 
curvature of the ground state^ii. 



'-'T^q^q'- 



A3 



-P 



2 2 d Eg (eg) 



(50) 



This property is a direct consequence of the relation be- 
tween current components (|47p 



diz 



= 21 



X 



Xi 



which seems to be peculiar for the flux qubit only. For 
other types of superconducting qubits there is no simple 
relation between the tank detuning and the curvature of 
the ground state. 

As is seen from ([^^ , our theory predicts the modifica- 
tion of the tank quality factor due to the interaction with 
the qubit. In principle, this effect allows us to measure 
the dephasing rate T^p. However, for the parameters of 
the flux qubit used in Ref. this effect was rather weak 
to be experimentally detected. 



F. Low frequency response of the charge qubit 

For a charge qubit with different critical currents, Ici 
and 7(72 of its Josephson junctions the energy gap is 



(51) 



where Cj = 



Ej2 



2EjiEj2Cos4>x^ Eji = 
$o/ci/2^, Ej2 ^'$oW2^, C = AEc{l-ng). Here 
the polarization charge on the island Ug is controlled 
by the gate voltage Vg via the capacitance Cg, namely 
"9 = CgVg/e. 



The components of the current operator are as 
follows^: 



LX 



Iy = 



2ii EjiEj2 . , 

iz = : smfflv-, 

^ $0 2AE 



2ttEji+Ej2 . 0x 

^ sm 

$0 4e2 2 

/r, 7^ ^2 4:CEjiEj2 2 

[Eji - Ej2 ) T cos 



2tt Eji - Ej2 4>x 

^ r^5 cos — — 

$0 4e2 2 

t-17 I 17 ^2 , '^CEjiEj2 . 2 0x 
(Eji+Ej2) H : sm 



(52) 



(53) 



(54) 



In case of symmetrical junctions {Eji — Ej2 — Ej) we 
have: Iy — 0, 



where 



^^ = ¥;;a^^^"'^^' 

TiEjC . 4>x 

ix = r-;^ sm , 

$oAS 2 ' 



A£;= a/4£;2cos2I^ + C2. 



(55) 
(56) 

(57) 



For the tank detuning and friction we obtain from (|43p 
and (|44ll . neglecting the terms which are proportional to 
dp'^'^iyd^x- 

-J n L rtUJrp 



- UJ 



AE 



4) 



(58) 



7t 



4PLghWj.r^ (59) 



(A£;)3 

As distinct from the flux qubit, the tank detuning caused 
by the charge qubit cannot, in general, be expressed in 
terms of the ground state curvature. It can be done only 
for the symmetrical case (Eji — Ej2) at the point <j)x — 



From expression (|58p we calculate the dependence of 
the phase shift a of the output signal on the gate charge 
parameter Ug for different values of magnetic flux <^x /'^o 
applied to the qubit loop (see Fig.[T]). For the calculations 
we take the following values for the tank, Qt = 1000, 
k'^ = 0.01, LOT = 27T X 50 MHz, and for the charge 
qubit, Lq = 1.5 nH, Eji/h = 25 GHz, Ej2/h = 29 GHz, 
Ec/h = 3.5 GHz. 

In conclusion to this section it is worth noting that 
the average current in the qubit loop (Iq) is proportional 
to the difference of equilibrium populations between the 
two qubit states only in the absence of interaction be- 
tween the qubit and the tank circuit: (iq) = —Izp^'^'^^- 
The other components of the current, Ix and Iy which 
give the main contribution to the tank response (the last 
terms in (I43p and (|44p '). appear in the average current 
only due to the interaction of the qubit with the tank. 
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FIG. 1: The dependence of the phase shift a on the gate 
parameter rig for different values of magnetic flux $x / <E>o ap- 
phed to the qubit loop. 



III. CHARACTERIZATION OF IRRADIATED 
QUBIT BY RABI SPECTROSCOPY 



The essence of the Rabi spectroscopy method is the 
following. Under microwave irradiation, which frequency 
is close to the gap frequency, the level structure of the 
global system (qubit, radiation field and their interac- 
tion) is a ladder of pairs of states where the spacing 
between two levels in the pair is equal to the Rabi fre- 
quency, the value of which is controlled by the power of 
microwaves^. Normally, Rabi oscillations are damped 
out with a rate, which depends on how strongly the sys- 
tem is coupled to the environment. However, if a sec- 
ond low-frequency source (in our case this source is a 
low frequency excitation from the tank) is applied si- 
multaneously to the qubit it responds with persistent 
low frequency oscillations. The amplitude of these low- 
frequency oscillations has a resonance at the Rabi fre- 
quency and its width is dependent on the damping rates 
of the system^^'^^. Therefore, the interaction of the qubit 
with the tank induces the transitions between the Rabi 
levels which result in low frequency oscillations of the 
current in a qubit loop which, in turn, result in the mea- 
surable response of the tank^^. In particular, the signa- 
ture of a high frequency resonance can be read out from 
the low frequency voltage-flux characteristic of the tank 
circuit. 

Below we assume that the interaction between the 
qubit and microwave field does not influence the dephas- 
ing and the relaxation rates. As was shown in Refs. [sol 
andlsil. this assumption is valid for relative weak driving, 
sufficient short correlation time of the heat bath Tc, and 
in the large temperature limit: F << A^, H/tc, ksT—. 



A. Low frequency linear susceptibility of an 
irradiated qubit coupled to a tank circuit 

In this subsection we summarize in a concise form the 
main results from Ref. [29l which are relevant for specific 
qubit applications. 

The Hamiltonian for an irradiated qubit which inter- 
acts with a low frequency tank circuit is as follows: 

Aii' 

H = —az + hwoa+a + ff^t + H\l^ , (60) 

where the second term in ([SO]) describes a microwave field, 
the third term describes the interaction of the qubit with 
this microwave field: 

H^m^-\%F{a+ +a), (61) 

where Iq is the current operator of the qubit as given 
in ([22|) . F is the amplitude of the microwave field with 
the dimension of a magnetic fiux. The last term in ([50]) 
describes the interaction of the qubit with a tank circuit: 

H\l^ = Mlrlq = MIt {Izaz + Ixax + lyoy) • (62) 

We denote the eigenf unctions of the photon field as 
|iV): a+|7V) = ^JW+\\N +\), and a\N) = VN\N - 
1). The eigenf unctions of the noninteracting qubit and 
associated photon system we take in the form of a tensor 
product \1,N) = \1)(E)\N), and \2,N) = |2) ® \N). 

If the photon frequency luq is close to the qubit gap 
frequency AE/h, and the high frequency detuning S is 
small, 6 = ujQ — AE/h <C ujQ,AE/h, where for defini- 
tiveness we assume ^ > 0, then the energies of the 
states \1,N + 1) and |2,iV) are close to each other: 
Ei^N+i ~ E2,N — fiS. The same is true for the pairs 
of states \1,N) and |2,7V- 1), |l,iV + 2) and |2, -I- 1), 
and so on. Therefore, the energy levels of the system un- 
der consideration are a ladder of pairs of states which are 
specified by the photon number N. The spacing between 
two levels in the pair is equal to HS, and the distance 
between neighboring pairs is equal to the photon energy 
hcufr^. These levels of uncoupled qubit-photon system 
are modified due to the interaction This interac- 

tion results in an increase of the energy gap between two 
levels in the pair. The spacing between the states \1,N) 
and |2, — 1) becomes equal to hilu, where fij? is the 
Rabi frequency^^ 

nn = ^6^ + f)2, (63) 

where SSli = F^Jl\ + I^. 

From now on we will call these two nearby dressed 
states Rabi levels. 

As was shown in Ref. |29| the interaction ((62|) between 
the qubit and the tank results in the transitions between 
Rabi levels. It is just these transitions which result in the 
low frequency response of a qubit detected by the tank. 
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In Ref. |29| we obtained the evolution equations for the 
elements of the density matrix which describe the tran- 
sitions between these Rabi levels: p, p+, and p_, where 
p is the difference of the populations between higher and 
lower Rabi levels. (Note, that here the definition of p is 
different from that given in Section II) . These elements of 
the density matrix are usually accounted for by a so called 
rotating wave approximation (RWA). The equations for 
the elements of the density matrix which describe in RWA 
the interaction of an irradiated qubit with the tank are 
as follows: 



dp 



-Aip+Bp^ 



2Trh 



Iz sm20p^+T^cos20, (64) 



-Bp-A2p+- 



2'Kh 



B 



1 



sin 29 cos 29. 



(69) 



The angle 9 is given by tan 29 = 
TT, so that cos 20 ~ cos I 



smf 



72 + iL) 



1/2 



-Qi/S, where < 26* < 



The steady state solution for the elements of the den- 
sity matrix and the low frequency linear susceptibility of 
a qubit are as follows2^: 



/zC0s2e'p_+r_ sin 29, 
(65) 



JO) 



(r2 



Y2 

Ti 



r_ cos 26*, 



(70) 



= -zr^flP+-ry/9_- ^^^^° /z (p+ COS26' - psin26') , 



(66) 



p2 



r_sin26l, 



(71) 



A. 



Ti 



cos^ 29 -t- sin^ 29 



— sin^ 29 + cos^ 261 



(67) 



(68) 



(0) .^R (0) 



(72) 



, s 2$or2_R J (0) 



sin 26* 



02 

—^B cos 261 



(73) 



2$of^i? 



^^+(")^' i?("J)2;ir/ ^^+^^°^'^ 



1 \ 

(zw + r^) ( + ^ ) - + ^i) Y' 



(74) 



2$o^^ 



D(uj)2TThTl 



IzpX cos 20 ( + — ) (iw + 2r^) , 



(75) 



where 



D{u) = {ioj + T^Y ( + ^ ) + (icj + Ai) 



(76) 



It is interesting to note that under high frequency ir- 
radiation the population of the Rabi levels becomes in- 



verted. It is seen from (iTOl) . where the quantity p' 



which is, by definition, the difference of the populations 
between higher and lower Rabi levels, is positive, since 
for (5 > we have cos 26* = —6/Q.b. < 0, and always 



r_ < 0. In addition, as S tends to zero, which 

causes the equalization of the population of the two lev- 
els {pii — P22 = 5) when the high frequency irradiation 
is in exact resonance with the energy gap of the qubit. 
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B. Calculation of Zi(tj) and Ziiyi) for irradiated 
qubit 

Since the matrix elements of ax and cry between the 
Rabi levels are zero, the average current is: 



(Iq) = IzTr{Gaz) ^ Iz{p{t) cos29 + p+{t) sm29) . 



By using Eqs. - (p^l) we obtain: 
at ±1 

+ mRTiS\n2ep^{t)]+ IzT^ 
This equation can be rewritten as follows: 



(77) 



1 



= - iIznRsin20p^{t) + IzT^ 

at 1 1 



(78) 



(79) 



With the aid of the steady state solutions ([70]) - (|7^ we 
find the stationary current when the interaction between 
the qubit and the tank is absent: 



1 + 



1 + 52T| + TiT2nl ' 



(80) 



where T2 = l/T^. 

It is not unexpected that this expression is quite similar 
to the one for the longitudinal magnetization in NMR^i. 
However, the important difference is that the frequency 
ill depends essentially on the current components of the 
qubit {rii =F^ll +1'^). 

For the Fourier component of (|78pwe obtain: 




27rAf , , Jz r / ^ 

= -^^T{uj)7f- [XpHcos29 

q?o J-i 

+ Xp+ {^) sin 29 + sin 29xp_ (w)] ■ (81) 

Therefore, for irradiated qubit the quantities Zi{lu) and 
Z2{u>) are as follows: 



Zi{lo) 



[x;(w)cos20 



Xp+ {^) sin 29 - nRT,Xp_ (lo) sin 29 



(82) 



Z2{uj) = -^^^[x;{co)cos29 

+ Xp+ {^) sin 29 + nRT,x'p_ (lo) sin 29] . (83) 

In this case the tank detuning ^ and the friction Tt can 
readily be obtained from (|4]) and (O by using the expres- 
sions ([83l) and ((82|) for Z2{uj) and Zi{uj), respectively. 



91,2 r r2, ,2 



ujhnRTi T^ + s^ + TiT^nl 



-fp (^) + j^/p+ (^) /p- (^) 



(84) 



IT- P 



9L2r r2 2 



;2ftr!flTi r2 + ,52 + r^r^f|2 



"/p(^) + 7^/p+ (^) + — r ^p- 



(85) 



where 



D{oj) 



sin 26* 



(iw + Fy) ( iw-l- ^ 



Q2 

-^B COS 29 

F„ 



171 



(86) 



(iw + Fy) ( ^ 



f]2 

(ic^ + Ai)-^ 
i if 

(87) 



— ) (^w + 2F^) . 



/p-H 



The expressions ((84)) and (|85|) can equally be used for 
the flux and the charge qubit with due account for the 
different current components ( (|47p for the flux qubit and 
([5^ - ([M)) for the charge qubit), and the different struc- 
ture of the energy gap f (|45|) for the flux qubit and ((5T|) 
for the charge qubit). 



C. Damping and amplification by a 
superconducting flux qubit 

A remarkable property of the irradiated flux qubit cou- 
pled to a low frequency LC circuit is that the effective 
quality factor of the tank can be made several times 
higher (lower) than that for an unloaded tank. A clear 
phenomenological explanation of this effect, which has 
been found experimentally, has been given in Ref. [25l . 
where the effect has also been explained by a rigorous 
quantum treatment of the problem followed by a numer- 
ical solution of the corresponding equations for the global 
density matrix of the qubit-tank system. 

Here we show analytically that in the frame of our ap- 
proach the effect results from the frequency dependent 
quality factor of the tank coupled to an irradiated qubit. 
The analytical expression we obtain for the quality fac- 
tor is valid for small high frequency detuning, S << wq, 
AE/h which allows us to account only for the transi- 
tions between the Rabi levels. As is seen from (1551) the 
amplification occurs if the high frequency detuning 5 is 
positive. It means that on the energy diagram of the 
qubit the working point of the driving current from the 
tank is located to the left from the point of the high fre- 
quency resonance (see Fig. 1 in Ref. [25h . Positive values 
of S correspond to the inverted population of the Rabi 
levels (see Eq. ([70)) ). It is just the emission of these Rabi 
photons which pumps the energy into the tank. For neg- 
ative values of S the working point of the driving current 
from the tank is located to the right from the point of the 
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FIG. 2: (Color online). The dependence of the quality factor 
of the tank coupled to an irradiated qubit. The quality of un- 
loaded tank Qt = 300. The upper curve { fx = 0.0027, 5 > 0) 
shows the effect of the amplification. The lower curve [fx = 
0.004, 5 < 0) demonstrates the effect of the damping. 



high frequency resonance. For this case the lower Rabi 
level is more populated than the higher level, and the 
tank transfers the energy to the Rabi levels. 

From (|55|) the quality factor of the tank coupled to the 
qubit can be written as i = 77 — \- 77—, where 



1 

Qqb 



52 + T^V^Vl\ 



(89) 



As is seen from Q — Qt at the point fx = Q {e = 
2^oIqfx)- The reason for this is the vanishing of the 
current component Iz at this point (see (|47p). 

We calculate from expression (1851) the frequency de- 
pendence of the tank quality factor Q and the voltage- 
frequency curve for the tank coupled to the irradiated 
qubit. We take the following parameters of the tank, 
Qt = 300, fc2 = 0.001, ljt = ^tt x 50 MHz, and of the 
qubit, Lg = 25 pH, /, = 300 nA, = 4 x 10^ c-\ 
Ti = 1.25 X 10-8 c, A/h = 3.78 GHz. For the microwave 
frequency we take wq = Stt x 3.81 GHz, and the mi- 
crowave power in the frequency units is F/h = 27r x 45 
MHz. 

The dependence of the tank quality factor on the low 
frequency is shown in Fig. [51 

The dependence of the normalized tank voltage on the 
low frequency (bias frequency of the tank) calculated for 
the same data as those for Fig. [5] is shown in Fig. [31 
For the same data we also show in Fig. [3| the depen- 
dence of the voltage and the tank quality on the bias 
DG flux fx- The curves are plotted for a tank resonance 
uj = LOT- The blue (short) arrows at this figure show two 
points where the microwave frequency matches exactly 
the gap frequency wq = A^/S. Between these points the 



FIG. 3: (Color online). The dependence of the normal- 
ized tank voltage on the low frequency. The quality fac- 
tor of the unloaded tank Qt = 300. The upper curve 
(fx = 0.0027, S > 0) shows the effect of the amplification. 
The lower curve {fx = 0.004, <5 < 0) is the effect of the damp- 
ing. The curve, which is in between, is the voltage curve for 
the unloaded tank. 
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FIG. 4: (Color online). The dependence of the normalized 
tank voltage (black curve) and the tank quality factor (red 
dash-dot curve) on the bias flux fx, plotted for uj — ujt- The 
quality factor of the unloaded tank {Qt = 300) is shown by 
the green straight line. The black (long) arrows show the 
points which correspond to the voltage at the tank resonance 
in Fig. [3] The blue (short) arrows show the points where 
the microwave frequency matches exactly the gap frequency 



high frequency detuning S is positive, and the Rabi pho- 
tons pump the energy into the tank. Beyond them the 
high frequency detuning S is negative, and the energy is 
drained off the tank. 

The analysis of Eqs. ([M)) and (|55|) shows that the de- 
tuning of the tank due to the interaction with the qubit 
is rather small, therefore, the resonance frequency of the 
coupled tank is practically equal to the unloaded value 
lut- However, the quality factor of the coupled tank has 
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FIG. 5: Flux qubit coupled to the inductive DC SQUID, 
which is incorporated in a resonant tank circuit. 



a strong dependence on the frequency. Generally, the 
position of the peak value of the quality factor does not 
coincide with the resonance frequency of the tank (see 
Fig. [2]). Therefore, the tank voltage curve can be appre- 
ciably differed from the lorentzian line. 



IV. CONTINUOUS READOUT OF A 
SUPERCONDUCTING QUBIT USING A DC 
SQUID IN THE INDUCTIVE MODE 

This method has been realized in Refs. [l^ and [H. 
In these works a DC SQUID operating in the inductive 
mode (the bias current through the SQUID is less than its 
critical current) was inductively coupled to a flux qubit. 
To enhance the measurement sensitivity, DC SQUID was 
incorporated in a resonant LC circuit (see Fig . [S]). The 
spectroscopic measurements in Refs. fisl andlssThave been 
performed by measuring the switching currents of the DC 
SQUID after the microwave pulse had been applied to the 
qubit. 

Here we propose the modification of the method which 
does not require a measurement of the switching cur- 
rents. Our method is based on the detection of the re- 
sponse of the resonant LC circuit in which a DC SQUID 
is incorporated. It is assumed that in the process of the 
measurement the qubit is irradiated by microwaves con- 
tinuously. We show that the DC SQUID in the inductive 
mode directly measures the Iz component of the current 
in a qubit loop. Our approach can be applied both for 
flux and charge qubits. 

First, we write the well known expression for the cur- 
rent across the DC SQUID 



II = /cC0s(/5,)sin7, 



(90) 



where 1^ is the current in the inductor Lt, which is con- 
nected in series with the SQUID loop, Ic is the critical 
current of the SQUID, 7 is the Josephson phase, which is 
related to the voltage across a SQUID by the Josephson 



expression Vsq = fj-^- The quantity fgq is the normal- 
ized total flux in the SQUID loop, $sq, which is the sum 
of the external control flux and the flux M{Iq) from 
the qubit: 



$0 



7r$ 



X 



$0 



7rA/(/g) 
$0 



(91) 



In this expression we neglected the flux Mldrc, which 
is generated in the SQUID loop by a circulating current 
Icirc- Below we assume II << Ic, which implies 7 << 1. 
From ([90)1 we get 7 = Il/Ic co^ fsq, and for the voltage 
across the tank we have: 



V=iLT + L, 



diL 
dt 



where 



2-kIc cos fsq 



(92) 



(93) 



is the Josephson inductance of the DC SQUID. The equa- 
tion for the voltage across the tank is similar to ([1]): 



V + -irV 



1 



-V 



ijjipLTlh- 



(94) 



We assume the coupling of the qubit to the tank is weak: 

'^^^j^^"^ << 1. For this case Eq. ([94]) can be written in 
the following form: 



V + 7tU -f ( 



cos/ 



X 



-V 



cos I 



X 



sin/. 



X 



cos/x + -r- 



-V = ujlLrh, 



(95) 



where fv - "''^^ - ^0 

wneie /X - , - 2-kIc ' 

As it is seen from this expression the tank response is 
proportional to the averaged qubit current, but not to its 
time derivative. Therefore, we may take for the quantity 
(iq) in ([55]) its stationary time independent part (iqY^- 
By doing this we may neglect the off diagonal components 
of the density matrix , and p- since their stationary 
parts are proportional to the coupling constant fc, and 
their contribution to the tank response scales as . 

Therefore, the influence of the qubit results in the tank 
detuning 



\Lt 



Ic ' 



(96) 



where we take fx = 7r/4 and assume i^^"'' << Lt- The 
detuning can be detected from the phase of the output 
signal at resonance: 



tana = — 



%/2 



Lt 



{I, 



Ic 



(97) 
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The expressions (|96)) and ([97]) can equally be applied both 
for ground state and spectroscopic measurements. For 
the ground state measurements (the microwaves are ab- 
sent) we have from (jSTj) (/g)^* = —Izp^'^'^^ For spectro- 
scopic measurements of an irradiated qubit (Iq)^* is given 
in dSOl). 

In fact, DC SQUID measures not an absolute value of 
the flux threading its loop, but the change of the flux. 
Therefore, we may subtract the ground state measure- 
ments from those for an irradiated qubit to obtain a pure 
lorentzian line 

f^- \ ft- \ ( Lq\^ Iz (en) 

(tana)^,-(tana)_ = -^(^ — j -p^ 

(98) 

which allows for the direct measurement from its width 
of the dephasing time T^^^ . 

This technique is similar to that used in low field NMR, 
where a longitudinal magnetization can be measured di- 
rectly with the aid of SQUID (see Ref. IIO and references 
therein) . 

In conclusion, we developed here a quantum theory for 



the calculation of the low frequency linear susceptibility 
of the interferometer type superconducting qubits. The 
obtained general results are applied for the calculation of 
the tank detuning and friction both for flux and charge 
qubit. For irradiated flux qubit wc obtained explicit ex- 
pression for the tank quality factor which allows us to cal- 
culate the recently found effect of the amplification and 
the damping of the tank. We have also shown the ap- 
plication of the theory to the continuous radio-frequency 
monitoring of the qubit with the aid of DC SQUID in the 
inductive mode. Our theory shows that radio-frequency 
method can also be applied for the investigation of other 
types of quantum two level structures provided the inter- 
action between the tank and quantum system is weak. 
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